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Abstract 

The phase diagram of the model of spinless fermions with repulsive nearest neigh- 
bour interaction is calculated analytically on a hypercubic lattice in infinite dimen- 
sions (d oo). In spite of its simplicity the model displays a rich phase diagram 
depending on the doping 5, the interaction U and the temperature T. The system 
can be in the homogeneous phase (hom), the nonsegregated AB charge density wave 
(ab-CDW), the AB phase separation region (ps-AB/hom; coexistence of AB-CDW and 
hom), the incommensurate phase (ip) or the IP phase separation region (pS-AB/iP; 
coexistence of AB-CDW and IP). We identify three important values of the interaction 
t^iPL = 0.572 < Uipn = 1.914 < Uip/ps = 4.212 which distinguish four intervals of U. 
These imply four different types of phase diagrams. In all the three phase diagrams with 
U below Uip/ps the IP appears. We propose a new general ansatz for the order parameter 
of this phase. A competition between the IP, the PS-Ab/ip and the PS-Ab/hoM is found. 
The relevance of our findings for the phase scenario of the Hubbard model is shown. 
^Present address: Laboratoire de Physique des Solidcs, Universite Paris-Sud, bat. 510, Centre Univer sit aire, 
F-91405 ORSAY Cedex 
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1 Introduction 



The interest in non-perturbative results for highly correlated itinerant fermion systems revived 
during the last years. Although in general analytically exact solutions of fermionic models 
with strong interaction are beyond reach, particular cases in the parameter space are solvable. 
The dimension plays an important role in this respect. In dimension d = 1, for example, the 
Bethe ansatz is at our disposal by which numerous models are analytically tractable. 

The other extreme, i.e. d— )>oo, is a useful tool for classical systems and has been extended 
to itinerant quantum systems about six years ago [1, 2]. This technique is successful for a 
variety of itinerant quantum models such as the Hubbard model [3-8], the t-J model [9, 10], 
the Anderson model of disorder [11], the periodic Anderson model [12] and the Falicov-Kimball 
model [13, 14]. An extensive review is written by Vollhardt [15]. 

In the present paper we study the model of spinless fermions with nearest neighbour 
interaction on a liypercubic lattice in the limit of infinite dimensions. It is our goal to find 
the phase diagram in the ground state as well as for finite temperatures. Our analysis shows 
a competition between the homogeneous phase (hom), the nonsegregated AB charge density 
wave (ab-CDw), the AB phase separation region (ps-ab/hom), the incommensurate phase 
(ip) and the IP phase separation region (ps-ab/ip). Phase separation means the segregation 
of phases: one part of the system has a doping above the average, another has a doping below 
the average. In the case PS-AB/hom, a hole-rich homogeneous phase and a hole-poor AB 
charge density wave coexist. In the case PS-AB/ip, a hole-rich incommensurate phase and a 
hole-poor AB charge density wave coexist. These phenomena are for the first time comprised 
in a single model simple enough to be analysed analytically. 

Some of the phases mentioned above attract attention since the discovery of the high-T^ 
superconducting materials. Phase separation is a manifestation of an attractive tendency. 
This is found for example in the t-J model [16-19] and in effective Heisenberg models [20]. 
Experimentally this is observed in lanthanum cuprates [21]. It is suggested that bound states 
of two or more particles will occur close to the region where the phases separate [16, 17, 19]. 
Spontaneous symmetry breaking with incommensurate order parameter was found theoret- 
ically in itinerant quantum antiferromagnets [22] for the Hubbard model [23-27], the t-J 
model [28-30], and the Falicov-Kimball model [31]. Experimentally it was observed that the 
incommensurate order parameter couples to the superconducting order parameter [32]. 

The model of interacting spinless fermions can serve as a relatively simple model which 
allows to test certain approaches like the limit of infinite dimensions. Moreover it allows to 
obtain an overview over the possible phase scenarios in more elaborate models and the com- 
petition between these scenarios. As we will see in the course of this paper many qualitative 
features of the infinite dimensional spinless fermion model are encountered in the infinite 
dimensional Hubbard model and other models, too. 
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Furthermore there are physical systems which can be described with the model of interact- 
ing spinless fermions such as strongly polarised systems, e.g. a lattice gas description of a ^He 
in a strong magnetic field, or ferro (ferri) magnetic electronic systems where one spin-band is 
filled. In these systems only electrons of one spin species are itinerant. The latter situation 
is realised in magnetite (Fe304), where the lowest spin- up band is half filled, leading to a 
metallic conductivity above a temperature ~ 119K at atmospheric pressure. At the 
system undergoes the Verwey-transition into an insulator [33]. Cullen and Callen [34, 35] first 
suggested a model of spinless fermions with nearest neighbour interaction to describe this 
transition. 

The model itself can be solved exactly in d = 1 by Bethe ansatz [36, 37]. At half filling 
the solution describes a transition from a homogeneous phase to a AB charge density ordered 
phase at a finite value of the interaction. The transition is discussed in detail by Shankar [38], 
who also calculates various response functions. A transition to a charge density ordered state 
at a finite interaction turns out to be a peculiarity of one dimension [39]. For d > 1 at half 
filling on a bipartite lattice the transition is shifted to arbitrarily small interactions according 
to renormalisation group results [39, 40]. Using a quantum Monte Carlo simulation technique 
the problem was studied in two dimensions by Gubernatis et al. [41]. For finite temperatures 
and half filling in d > 1 Lemberger and Macris [42] proved rigorously the presence of an 
AB-CDW at sufficiently large values of the interaction. Apart from Yang and Yang [37], who 
include certain aspects of finite doping, the research so far concentrates on the half filled case. 

It is the objective of the present paper to provide an exhaustive description of the phase 
diagram of repulsively interacting spinless fermions thereby extending and compiling previous 
results [43, 44]. The striking analogy between certain results concerning the incommensurate 
phases in the model of spinless fermions and in the extended Hubbard model [45] triggers the 
interest in a circumstantial discussion of such calculations. 

The paper is organised as follows. In sect. 2 the Hamiltonian and the relevant results for 
d — > oo are given. Subsequently in sect. 3, the equations are derived which determine the 
second and the first order transition lines at zero temperature. From these the T = phase 
diagram is constructed. The phase diagram for finite temperatures is constructed in sect. 4. 
In sect. 5 we conclude our work by a summary and a discussion of the implications of our 
results for the understanding of the phase scenario in the (extended) Hubbard model. 

2 Hamiltonian and limit of infinite dimensions 

We study the model of itinerant spinless fermions with nearest neighbour hopping t and 
nearest neighbour repulsion [/ on a hypercubic lattice in d^oo dimensions. The Hamiltonian 
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(1) 



where cl (cj are the creation (destruction) operators for ferniions on site := cfcj and 
U--, t-- are zero if i,j are not nearest neighbours. The hopping matrix elements are scaled 



and j are nearest neighbours. This scaling guarantees that the model remains non-trivial in 
the limit d^oc [1, 2]. Then the dispersion relation is 



Since the hypercubic lattice is bipartite, the Hamiltonian is particle-hole symmetric and we 
may restrict the filling n to < n < 1/2 or, equivalently, the doping 5 to S := 1/2 — n > 0. 
Throughout the paper we will work in units oi h = 1, = 1 and t = 1. The lattice spacing 
is set to unity. 

It is known that in the limit of large dimensions the self-energy becomes site-diagonal and 
thus k independent [15]. For interactions as the one considered in (1) which are not on-site the 
Hartree approximation becomes exact [2]. We show that on the two-particle level the random 
phase approximation (rpa) becomes exact on d— > oo. The proofs of these statements rely 
all on the resummation of the corresponding perturbation series on a diagrammatic level. In 
order to cover both the Hartree approximation for the one-particle quantities and the random 
phase approximation for the two-particle quantities we argue on the basis of the generating 
functional $[G] introduced by Baym and KadanofF [46, 47]. Furthermore it is advantageous 
that this line of argument can easily be extended to 1/d corrections [48]. 

The perturbation diagrams of which the sum equals the generating functional $[G] are 
classified according to their contribution in powers of 1/d. In fig. 1 three examples are given. 
The diagram in fig. 1(a) is of order d^ whereas the diagrams in 1(b) and 1(c) are or order d~^ . 
The power counting for the diagrams relies on the proportionality of the propagator in real 
space G- oc [1^ 2] where the Manhattan metric || • || counts the minimum number of 

hops to reach fj starting from f-. The proportionality can easily be derived by expanding the 
propagator in the scaled hopping. For instance, in the diagram in fig. 1(a) both propagators 
are local (i = j). The factor 1/d from the scaling of the interaction is compensated by the 2d 
realisations of this diagram on the lattice (there are 2d nearest neighbours). In the diagram 
in fig. 1(b) both propagators go from one site to one of its nearest neighbours (||^i — ''jll = 1) 
which adds an extra factor 1 /d. 

The one- and the two-particle irreducible vertices are found from $[G] by functional 
derivation with respect to the Green function of the interacting system S = S^/5G and 
S = S^^/SGSG' . Inserting S (S) into the Dyson equation (the Bethe-Salpeter equation) 
yields the full Green function (the general two-particle correlation function). The ensuing 



as t-j := - 



t/\/2d and the interaction matrix elements as [/• = U/(2d), given that the sites i 
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diagrammatic series are just the self-consistent Hartree and random phase approximations. 
They are depicted in figs. 2(a) and 2(b). 

Two further consequences of d— )-oo are important in our calculations. First, the DOS of 
the non-interacting system is a Gaussian 

P(^) ■■= ^?'^(^-^a;) = -^^^p(-y) ' (3) 



k 

where N is the system size. Eq. (3) follows from the central limit theorem since the dispersion 
in (2) is a sum of terms each depending only on one component of k [1]. Second, the two- 
energy density p(ei, 63, k) which occurs in the calculation of response functions (see below) is 
a Gaussian, too. It depends on k only via the parameter rj^ [2, 49] 
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^^P 2^ with 



1 

"Ik ■= ^Ecos(A;J . (4b) 

n = l 

In all dimensions — 1 < < 1 holds. We will call the parameter of incommensurability 
since the commensurate values of the wave vector k which fulfil = r]g for — 1 < r]g < 1 are 
of measure zero in the Brillouin zone. For d— )■ 00 is zero for most possible values of k as can 
be seen comparing the scaling of with the scaling of in (2). But for certain values of k 
which are of measure zero in the Brillouin zone for d^oc rir does not vanish. These values are 

' k 

physically very important since the value of rj is fixed externally; it is impressed on the system 
from the outside. For instance, one has r]^ = —1 for the nesting vector Q := (tt, tt, . . . ,7r)^. 
In contrast to this, the dispersion e- is involved in integrated quantities. This explains why 
both quantities rj^ and are physically relevant in spite of their different scaling in 1/d. 



3 Phase diagram at zero temperature 

Two regions of the phase diagram are known. At half filling and large interaction the AB-CDW 
is present and for very low values of the interaction and large doping the system is in the 
homogeneous phase. In the first case, the interaction requires that there are no two particles 
adjacent to each other. At half filling this requirement can only be met in the ab-CDW. In 
the latter limit, the interaction becomes finally irrelevant. So the system is in the HOM phase. 

Starting in the HOM phase for low values of the interaction and for low filling one ap- 
proaches the transition boundaries to other phases by continuously increasing the interaction 
and the filling. At the boundary in the U — S plane at which a second order transition from 
the HOM phase into another phase occurs the system becomes infinitely susceptible to density 
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fluctuations characterised by a certain wave vector k. The divergence of the density-density 
response function x(A;) indicates the instabihty of the HOM phase. This stabihty analysis is 
carried out in the following flrst subsection. The resulting second order boundary marks the 
maximum region the HOM phase may occupy in the U — S plane. 

Second order transitions can be superseded by first order transitions. The boundary which 
marks a first order transition between two phases can be constructed only if the ground state 
energies of both phases are known. In the second subsection we will construct the first order 
transition line between HOM and AB-CDW as well as the the first order transition line between 
IP and AB-CDW 



3.1 Second order transitions 

In the HOM phase the Green function G and the density-density response function x ^re 
diagonal in fc-space. The Green function as function of the Matsubara frequency uj^ and the 
wave vector k is 

G{iuj^,k) = {iuj^ - ej: + jiyK (5) 
The only trace of the interaction is the renormalised chemical potential 

jj = fl — Save = I^ — IT-U . (6) 

It equals the true chemical potential reduced by the space-averaged self-energy which is nU . 
Integrating (3) yields the relation between jl and doping 8 

^ = -5-KA)- 

The RPA series for x(A;) derived from fig. 2(b) is shown in fig. 2(c). It is a geometric series 
with the factor —Urj^Xoi^)- This follows from the evaluation rules for two-particle diagrams: 
the wavy interaction line with the wave vector k stands for the negative Fourier transform of 
U^j which is —Ur]^. A single bubble represents Xo(A;), i.e. one has 

xCk) = H^^^^- (8) 
1 + Urj^^xoik) 



The value of Xo(A;) is explicitly given by 
T 



Xoik) ■■= -^Y.^„r,Mi^n,k' + k)G{iu;^,k') =^ (9a) 



— oo — oo ^ 

J J 6i €n 



-QC 



/ , e.J^t] . (9d) 



v/27r(l - v) { tp-{l+r,)t 
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Partial fraction expansion leads to (9b) and the two-energy density (4) is used to transform 
the fc-sum over the Brillouin zone. The Fermi functions /F(e) := (1 + exp(P(x — fl))^^ in (9c) 
are step functions at T = but at finite temperatures the formula holds in the form (9c). 

The central result of (8) and (9) is that % depends on k only via the parameter rj^ [43]. 

The second order boundaries are found from 

= xiv^r' = xo{v^)-' + Uv-, . (10) 

The case 77 = —1 corresponds to k = Q. Its divergence indicates the transition to the AB-CDW. 
The integral (9d) simplifies for rj = —1 and one has 

^ =Ei(-/iV2) (11) 



U 2v^ 

for the corresponding values of the interaction. Henceforth, we will drop the explicit notation 
of the /c-dependence of rj. 

The analysis of the relation (10) for general values of rj is similar to the one of Lindner and 
Schumacher [23] for the two dimensional Hubbard model and to the one of Freericks [31] for 
the infinite dimensional Falicov-Kimball model. Coming from low interaction and low filling 
we look for the minimum value of the interaction at which the susceptibility diverges for some 
value of 7]. We call the value of the parameter of incommensurability at which this happens 
r]„pt since it is the most favourable (optimum) value of rj for a phase transition. 

The results are depicted in fig. 3. The dotted curve in fig. 3(a) displays the line at which 
x(— 1) diverges. It is found from (7) and (11). The solid curve in fig. 3(a) displays the line 
at which x(r]opt) diverges. It is found from (7), (10) and the derivative of (10) with respect to 
7]. For the latter curve the dependence of ?y„p, on the doping is given in fig. 3(b). We observe 
that r]„pt((5) is continuous for 5 — > but jumps to —1 at 5 = 5end where the two curves in fig. 
3(a) intersect. The essential point is that for < 5 < S^^^ the system is more susceptible 
towards an incommensurate CDW than towards the usual AB-CDW [43]. 

The optimum value rj„pt does not exceed 7]„pt < —0.793 at zero temperature. This is directly 
related to the nature of the interaction. Since the repulsive interaction takes place between 
particles on nearest neighbour sites an AB-structure is favoured even by the interaction alone. 
On the contrary, a local interaction alone does not favour any specific spatial pattern. There- 
fore one may expect higher values of rj^^t in models with local interaction. Indeed, Freericks 
[31] finds values of rj^^t up to one, i.e. — 1 < rj^^t ^ 1, for the Falicov-Kimball model. For 
the Hubbard model in d = 3 Lindner and Schumacher [23] find — 1 < rj^^t < 0, just as Doll, 
Dzierzawa, Fresard and Wolfle [25-27] do. On the other hand, there are indications that the 
quantum fiuctuations in the Hubbard model prevent rj^^t from deviating too far from —1 [50]. 

The vicinity of the perfect nesting singularity at 5 = is of particular interest. In the 
limit S ^ the optimum value of the parameter of incommensurability rj is related to the 
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doping by 

l + r]opt = ^S^ + 0{S^) where (12a) 
2sD(s) = 1 ^ s ~ 0.924139 . (12b) 

The function D{x) is the Dawson function defined by D{x) := —^^/^^exp(—x'^) eri{ix). With 
(12) the relation between 5 and U in the limit 5 — )■ or equivalently in the limit U ^ can 
be derived from (10) and (7) 

= -^exp(-V2^/U -C/2 + S) , (13a) 

1 °° 

5 = ln(s) - f dtexp{-{t - sf) \n{t^) ^ 0.24536 , (13b) 

2 Wit J 



where C is Eulers constant: C ~ 0.57722. The corresponding relation for fixed value rj = —1 
is derived from (11) and (7) 

5^^ = -^exp [-V2^/U - C/2) . (14) 

The fact that S-^-p > Sj^-^ at fixed U shows that the divergence of the susceptibility for an 
incommensurate value of the wave vector k occurs before the divergence for the nesting vector 

Q. 

Meanwhile van Dongen found that the relations (12), (13) and (14) are also relevant in 
the extended Hubbard model in the weak coupling limit [45]. The main difference to the 
situation in the model of spinless fermions consists in a renormalisation of the interaction due 
to fluctuations. 

Using (12) and expanding (4b) for small deviations of k from Q we find 

-^\k-Q\^^S . (15) 
Vd s 

The deviation between the optimum wave vector k and the nesting wave vector Q charac- 
terising the AB-CDW is proportional to the doping. Shraiman and Siggia found a similar 
proportionality for the incommensurate pitch of a spiral order in a quantum antiferromagnet 
[22]. 

3.2 First order transition HOM^ — ^AB-CDW 

In this subsection we derive that the AB-CDW is stable only at half filling since it is unstable 
towards phase separation at finite doping. Furthermore we discuss a promising ansatz for 
the order parameter in the IP. First, however, we have to introduce the description of the 

AB-CDW. 
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In the AB-CDW particle densities on the two sublattices of a bipartite hypercubic lattice 
are different. Let us define n to be the average particle density and h to be the difference 
between the particle density on the B-sublattice and n 

n ■= (^aeA> + KeB»/2 (16a) 
^ := (^fies)-'^- (16b) 

In Hartree theory the average self-energy is nU which is taken into account by the renormali- 
sation of the chemical potential (6). A finite value of h induces an alternating field A exp(«(5) 
{Q = ("^j • • • ) '^)^ is the nesting vector) which is A = hU on the sites of sublattice A and — A 
on the sites of sublattice B. In this way the translational symmetry is broken and the Green 
function and the susceptibility are no longer diagonal in /c-space. The mode at k couples to 
the mode at A; + Q. Since A; + (5 + (5 = A;uptoa reciprocal lattice vector it is sufficient 
to consider 2x2 matrices in the ab-CDW in contrast to the scalars in the HOM phase. The 
Green function becomes block diagonal in fc-space for finite h 



iuj„ + €- + /i 



{ico^ + jlf - el - A2 

From n = (T/N) G^i i(^'^n? ^) from b = (T/N) J]^ En G^i 2(^'^rn ^) obtains 



(17) 



n = J Mu;)p{V^:P^)e{u' - A') , (18a) 



-OO 

oo 



= j-j^=Uuj)p{V^J^^^)@{u;'-A^). (18b) 

— oo 

Here Q[x) is the Heaviside step function. These two equations determine jl and A (6, respec- 
tively) for externally given values for n and U . 

Evaluating the RPA series in fig. 2(c) with the propagator in (17) yields a 2 x 2 susceptibility 
X- The instability of the AB-CDW is manifest in det(x) < for all 5 > 0. At half filhng the 
AB-CDW is stable. This instability can be traced back to the non-convexity of the ground 
state energy which leads eventually to a first order transition and phase separation. 

The ground state energy _E'(n, A) can be foimd by integration from 9i?(n, A)/9ri| ^ = A* 
and ^-©(n, A)/9A|„ = with the help of (18) or using the relation between the generating 
functional $[G] and the free energy or the ground state energy, respectively [47], (see also 
the corresponding subsection for finite temperatures). The result is 



E(n, A; 7]) = — + -n^ + j duj sgn(w)Vt^2 + A2p(w) • (19) 

— oo 

The dashed-dotted curve in fig. 4 depicts a generic result for the ground state energy in 
the AB-CDW at [/ = 1.2. The dotted curve shows the continuation in the HOM phase for 
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S > 5ji^_AB- The system will segregate into a spatial region with a hole-rich HOM phase and 
into a spatial region with the AB-CDW without holes. Theoretically one has to pass to the 
convex hull of the dashed-dotted curve in fig. 4 which is a straight line from the point at half 
filling (5 = 0) and the point at 5 = (5ps_ab (Maxwell construction). The points (S-ps-ab,U) 
result from 

The functions K^{x) are modified Bessel functions. Since (5ps_AB is in the HOM phase (20a) 
results from (11). For (20b) we exploited the explicit integrability of (18) at half filling, i.e. 
for |/i| < |A|. Explicit integration in (19) leads to the condition for the straight line 

exp {^\kA^\- exp (-^^^\ = V2^ (^Sl,_,, - il,s-..S,s-..) • (21) 



4 \ 4 J ' \ 4 J "V 2 J V2 

Numerical evaluation of (20) and (21) leads to the first order transition line (solid curve) in 
fig. 5. For comparison the second order transition line (dotted curve) X()^(~^) plotted, 
too. This curve consists of the points {S^-„_j^^, U) which lie by construction to the left of the 
solid curve (cf. fig. 4). A point in the region PS-AB /hom stands for a phase separated system 
consisting of a hole-rich homogeneous phase (volume fraction v^) and a hole-free AB-CDW 
(volume fraction 1 — v-^) such that the average doping equals 8, i.e. = 5/5ps_ab 

3.3 First order transition IP^ — ^AB-CDW 

Although for most values of the interaction U the first order transition line turns out to lie to 
the right of the second order transition line for optimised parameter of incommensurability 
there is an intervall U G [f^iPL, f^ipn] = [0.572,1.914] where the situation is inverse (cf. fig. 
6). Approaching half filling from large doping xiVopt) diverges in this interval before the first 
order transition HOM^ — )■ AB-CDW occurs. The existence of this interval was not yet found in 
ref. [43] due to an incorrect evaluation of the Maxwell condition. 

The occurence of [/jpl can be understood already in the limit S 0. Expanding (21) to 
order S'^^_^^ using K-^^(x) = 1/x + 0{x) yields = 47r5pg_^g(l -|- C//V27r). From this one 
deduces with the help of (20) 

1 exp(-A/2^/C/-C/2 + ln(2)/2) 



^i + c//v^ 



(22) 



Comparing (22) with (14) and (13) for very small U one realizes that the following sequence 
holds: (5ps_4B > <5ip > 5ab- This implies that phase separation dominates for U — > at T = 
[43]. Exactly the same observation was also made by van Dongen in the extended Hubbard 
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model because the fluctuation induced renormalisation factor is equal in fact for all three 
relation (13), (14) and (22) [45]. 

The denominator in (22) leads to a change of the sequence at around U »i 0.56: S-^p > 
5ps-AB > <^AB- Thus one understands why in the exact calculation the incommensurate phase 
dominates the system above some particular value [/jpl ^ 0.572 of the interaction. Our 
calculations revealed that the sequence changes back to S^s-ab > <^ip at a higher value t/jpH ~ 
1.914 of the interaction and stays like that for U — )■ oo. 

If the system is in the incommensurate phase for U G [t^iPL? ^^iph] at some finite doping 
and in the ab-CDW at half filling the question remains where the transition takes place and 
which kind of transition it is. To give an answer to this question one needs the ground 
state energy of the incommensurate phase, i.e. one has to know of which structure the order 
parameter is. Note that this is a question without a straightforward answer even though the 
interaction induces only Hartree terms [51]. In the following we will discuss a new ansatz 
for the spatial structure of the order parameter which goes beyond the assumption of charge 
density waves characterised by well-defined wave vectors. In contrast to the case of spin 
ordering in the Hubbard model it is not possible to resort to a spiral ordering with sufficiently 
high symmetry which can be resolved by standard methods. 

For the spatial dependence of the order parameter b(r) := (n-) — n we use the product 
ansatz 

= b.iju^ir^^^) (23) 

where the functions u,; take the values ±1 with equal frequency over the whole lattice and r^*^ 
is the component i of f. So, the spatial dependence of 6(r) in each direction is independent 
from the dependence in the other directions. This form is most suitable for the description 
in the limit d^oo. Note that no periodicity of the tij(r*')) is required. 

In appendix A it is shown tliat the ground state energy depends on b{f) only via the 
amplitude 6q and the relative frequency h on the whole lattice that b(f) has the same sign on 
the adjacent sites. Unless h is strictly unity the average of b(f) over the whole lattice vanishes. 
Note that h can take all values between and 1 in the limit d— >oo; in finite dimensions the 
values of h are multiples of l/(2d) on hypercubic lattices. The amplitude 5q and the relative 
frequency h are the only parameters which are relevant for the influence of a spatial structure 
of the form (23) on the energy. 

In order to compute the energy functional for b(f) most easily we exploit the freedom of 
choosing the functions u-. Keeping h constant we set u-{r^'^^) := 1 for i < hd and u-{r^'^^) : = 
(— l)''^ ' for i > hd. This choice is a particular realization of the relative frequency h which 
restores periodicity and can be characterised by the wave vector 

= ( Q,0,...,0 ,5,7r,^...,7rjt . (24) 

hd {l-h)d 
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With the help of (24) one can define a relation between rj and h, namely r]Q^ = 2h — 1. 

Since the ground state energy E of the IP characterised by (23) is independent from how 
the relative frequency h is realized one can compute E for the CDW belonging to Q;,. In this 
way one obtains the ground state energy for a large class of order parameters most of which 
do not display any periodicity from a simple calculation involving only 2x2 matrices. This 
calculation, which can be found in appendix A, yields with A = —rjUb^ 



E{n,A,r^) = — -—-+ dv ) ' \ 



sgn{v) Vtt 



(25a) 



„2 



= - [dv ^ '""^ erf (P) . (25b) 



The actual energy is found from (25) by minimisation with respect to A and rj. The generic 
result for U = 1.2 G [C/ipl, t/jpn] is depicted in the solid curve in fig. 4. Indeed, the solid curve 
lies always below the dashed-dotted curve which belongs to the AB-CDW (rj = —1). But the 
minimisation with respect to rj does not remove the non-convexity. This means that there 
will still occur phase separation. For illustration the dopings 5div-AB (divergence of x(— 1)) 
and 5ps_AB (start of AB phase separation) are depicted which are relevant when only rj = —1 
is considered. The doping 5ji^_ip depicts where the susceptibility xiVopt) diverges. A second 
order transition to the incommensurate phase occurs which is stable from (5ji^_jp to 5ps_jp. 
At 5ps_ip the first order transition from the IP to the AB-CDW sets in, i.e. one constructs the 
convex hull to the solid line. The convex hull to the dashed-dotted line is no longer relevant. 
The rj values which are relevant in the IP are very close to those shown in fig. 3 (b). So we 
do not display them in a separate figure. 

Summarising the results of this section we conclude that for U G [f-'ipL, t/irn] the system is in 
the homogeneous phase at large doping; on decreasing doping it passes to the incommensurate 
phase by a second order transition and finally by a first order transition to the AB-CDW at 
half filling. The resulting T = phase diagram is plotted in fig. 6. For U < Uip-^ or U > C/jph 
no IP occurs and the system passes by a first order transition directly from the HOM phase to 
the AB-CDW at 5 = 0. 



4 Phase diagram at finite temperatures 

In principle, the arguments of the preceding section can be used to derive a phase diagram 
at T > 0, too. Yet there arise some new aspects. The energy as function of the doping is 
smoothed by thermal fiuctuations so that the regions of phase separation are shifted to higher 
values of the interaction and one finds also a second order transition HOM < — ^PS-AB/hom. 
Furthermore the incommensurate phase becomes more important at finite temperatures [43]. 
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4.1 Second order transitions 



First we investigate the stability of the homogeneous phase, which is stable at low interaction 
and low filling, on increasing interaction and filling. Starting from (9c) one can perform 
one integration after the substitutions 63 = U\J2{1 — rj^) + rje-^ or = U\j2[l — rj^) + 7762, 
respectively. With the following representation of the Dawson-function 

oc 

DM = (47r)(-^^2)vp j eY:g{-u^) j [u + u) du (26) 

— CO 

(vp for principal value) one obtains 



XM = 7 "^y^ D(^v^(l-r^)/(2(l+r^)) (/,(-.) - /,(.)) . (27) 



As at T = we compute the value of the parameter of incommensurability r] at which x{v) 
in (10) diverges first. The results for x(— 1) and for xiVopt) at T = 0.1 are presented in 
fig. 7(a). Comparison with the corresponding fig. 3(a) at zero temperature shows that at 
T > the intersections are shifted to larger interaction on the divergence line of xiVopt)- The 
corresponding values of rj^^^ at T = 0.1 are depicted in fig. 7(b). At ^start the value of rj^^t 
goes continuously to —1 whereas at 5end the value of r]„pt jumps discontinuously back to — 1. 
Note that finite temperatures lower the values of rj^^t] its maximum value is attained at T = 0. 

The interval of incommensurability [^start, ^end] becomes smaller on increasing tempera- 
ture. At Ti,,,, ^ 0.1826 the interval ceases to exist with (5start — <5end- The values of (5start 
and of 5end can be found in fig. 8. The dashed and the dotted curve enclose the area in the T-S 
plane where the divergence of x{v) occurs first for r] > —1. The dotted curve represents the 
dopings where 77„p, approaches —1 continuously and the dashed curve represents the dopings 
where r^„pt jumps to —1. 

In contrast to the situation at T = there exists at T > a region of dopings around half 
filling for which the ab-CDW is stable. Thus a complete stability analysis must investigate 
where the AB-CDW becomes unstable. One possible instability can be calculated from a 
generalised susceptibility which is given by the 2x2 density-density response matrix xiv)- 
One has to sum the geometric RPA series in fig. 2 taking into account the matrix form of the 
propagator in (17). This yields the extension of (10) to finite order parameter 

det{x-\v,^)) = det(xo\v) + Uv[l^^])=0. (28) 
The matrix in the r^-term accounts for the fact that Vr ^ = —Vr which is relevant for the 

k+Q 'k 

mode at k + Q. The matrix elements of 



Xoiv, A) 



(29) 



X^(r],A) Xo'iv,^) 

are derived in appendix B. The determinant appears in (28) in order to detect when xiVi A) 
diverges, i.e. when one of the eigen values of X~^{v^ A) becomes zero. The AB-CDW is stable 
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if and only if both eigen values are positive. Starting at half filling and increasing 5 we look 
for the value of rj at which x(v^ A) diverges first. By this procedure the dotted curve in fig. 
9 is obtained. 

Coming from 5 = in fig. 9 one realizes that the solid line ab precedes the dotted curve. 
This solid line marks the second order transition from the AB-CDW with r] = —1 to the 
incommensurate CDW with ?y > — 1 at finite amplitude of the CDWs. Note that the instability 
of the AB-CDW towards a deviation of 77 from -1 is not equivalent to its instability towards 
the superposition with a CDW with infinitesimal amplitude and different value of rj. The 
latter instability is signalled by the divergence of x(?7,A). The former is determined from 
the minimisation of the free energy given in (33) in the next subsection in the limit ?] — )■ — 1. 
Relying on our ansatz (23) we find that the instability of the AB-CDW towards a deviation of 
7] from -1 is the relevant instability since it precedes the divergence of x(r], A). This modifies 
our results presented previously [43]. 



4.2 First order transitions 

In order to calculate first order transitions one has to construct the convex hull to a non- 
convex free energy as function of the filling. The free energy functional for the AB-CDW at 
T > is given by 

F{n,A) = —--n'^ + i_m + Tjduj\n{l-Msgniu)VlJ^TA^))p{uj), (30) 

which can be constructed by integrating dF{n, A)/dn\A = and dF{n, A)/dA\n = or 
by using the relation between the grand canonical potential given by F[n, A) — //n and the 
generating functional $ in the Baym/Kadanoff formalism [47] 

F{n,A)-iin = - ($ + TTr (In(-G) - SG)) . (31) 

Since the cusp in the ground state energy at half filling (see fig. 4) is rounded at T > the 
Maxwell construction requires to determine two dopings: < The larger one, refers 
to the hole-rich HOM phase; the smaller one, 5^, to the hole-poor AB-CDW. Between S-^^ and 
S2 the convex hull is given by a double tangent with the slope //ave- So the conditions 

1/2-S, 

J dni2{n) = F{l/2-S„0)-F{l/2-S2,A) = {S2-S,)i2j,^^ (32a) 



1/2-s, 



dF dF 

— (1/2 - 0) = //AVE ^(1/2 - ^2, A) = //AVE (32b) 



have to be fulfilled for phase separation between the HOM phase and the AB-CDW to occur. 
They determine A, 5^ and S2 for given interaction and temperature. 

A generic result is shown in fig. 10 at T = 0.2. (At relatively high temperatures it is not 
necessary to account for the appearance of the incommensurate phase as we will see below.) 
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Below U ~ 5.0 no phase separation occurs. We observe a second order transition between 
the HOM phase and the AB-CDW. The point a in fig. 10 is the critical point where the second 
order transition line meets the two first order transition lines given by U{5i) and U{52). A 
point in the phase diagram 10 at a certain doping S in the region PS-AB/hom represents a 
phase mixture where the HOM phase and the AB-CDW coexists in spatially distinct domains. 
The AB-CDW occupies a region of volume fraction v-^ and the HOM phase a region of volume 
fraction 1 — v-^^ so that 5 = v^5^ + (1 — f i)52 holds. 

It is instructive to keep track of the loci of the critical points in the T-5 plane. They are 
given from (32) in the limit 62 — 6-^-^ 0. Since they lie all on the second order transition line 
HOM-( — )-AB-CDW their interaction co-ordinate is given hy U = X~^{v = The solid curve 
in fig. 8 depicts the critical points in the T-5 plane. 

From the origin to point a in fig. 8 the critical points of the AB phase separation lie 
within the intervals where the incommensurate phase is important. Let us define the tripel 
(<^ip/ps, 7ip/ps, t^ip/ps) = (0.2668,0.1825,4.212) by the co-ordinates of this point a. Then we 
can make the following statements: 

• Beyond point a phase separation dominates. For T > Tjp/ps (or for U > t/ip/ps) the 
HOM^ — ^-AB-CDW phase separation prevails, i.e. essentially no incommensurate phases 
occur. In this sense the phase diagram in fig. 10 is generic. 

• For T < Tip/ps (or for U < C/jp/pg) the incommensurate phase occurs. It will mask the 
critical point of the AB phase separation. The generic situation is shown in fig. 9 for 
T = 0.1<Tip/ps. 

In order to investigate phase separation including the incommensurate phase given by the 
ansatz (23) the free energy functional 

F{n,A;v) = - -n^ + iin + T J dio - U{lo)) p^{lo) with (33a) 

P,M = / , exp (- [u'lh + v'l{l - h)] /2) - e{u, v)) (33b) 

-00 27r^/i(l - h) 

for general r] must be known. It is derived in appendix A. The minimisation of the free energy 
in (33) with respect to A and to rj determines these values. In analogy to the case T = 
we construct the convex hull to F(n, A; rj) by the Maxwell construction. Like at T = we 
find a first order transition from the AB-CDW at low doping to the IP phase at higher doping. 
A region of phase separation PS-AB/ip occurs which is shown in fig. 9 between the curves 
ad and be. Point b is a new critical point where the second order line ab between IP and 
AB-CDW splits into the two first order lines. 

Comparing the case T > in fig. 9 with the case at T = in fig. 6 one finds that the 
region in which the pure IP phase is stable is considerably enlarged by finite temperatures. 
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This will become even more obvious by examination of the phase diagram at fixed interaction 
U. Up to now we discussed the phase diagram in the U — 5 plane at fixed temperature which 
was the best way to understand the physical principles which govern the interplay of the 
different phases. On the other hand, it is experimentally very interesting to investigate the 
system at given interaction because the temperature is more easily tuned than the value of 
the interaction. The theoretical interest in the phase diagram at given interaction is kindled 
by the appearance of a re-entrant behaviour (see below). 

We do not want to present here the whole line of argument for the derivation of the phase 
diagrams in the T — 5 plane which is analogous to the one for the U — S plane. We only state 
that for the T — S plane the two limits in which the character of the phase is known are the 
limit T — )■ oo and the limit T — > 0; 5 — )■ 0. In the former limit the phase is homogeneous; in 
the latter limit the phase is the AB-CDW. They serve as the starting points of the stability 
analysis. 

The results are shown in figs. ll(a)-(d). Depending on the value U of the interaction there 
are four generic cases: U G [C/jp/pg, oo], U G [C/jph, t/ip/ps], U G [C^ipl, t/jpn] and U G [0, C/jpl]. 

One recognizes the meaning of C/jp/pg which was defined analogously to Tjp/pg with the help 
of fig. 8. For values of the interaction U > C/jp/pg as in figure 11(a) the incommensurability is 
totally absent. For U < C/jp/pg an IP appears. 

Depending on the value of U one encounters different scenarios for T ^ 0. From fig. 6 
one sees that for U > C/jph the system does not display the IP at T = 0. So it appears only 
at finite temperatures as depicted in fig. 11(b). For [/jpH > U > C/ipl, however, fig. 6 tells us 
that the IP appears together with the ps-ab/ip down to T = as shown in fig. 11(c). For 
C/iPL > U > in fig. 11(d) there is again a phase separated region (ps-ab/hom) at T = 
similar to the case for C/ip/pg > U > C/jph (fig. 11(b)). 

Inspection of figs. ll(a)-(c) shows that re-entrant behaviour occurs since the phase tran- 
sition lines are not always functions of 5. (In fig. 11(c) the effect is hardly discernible 
in the plot.) For instance, in fig. 11(a) at 5 < 0.3 the system is in the HOM phase for 
temperatures < T < 0.25 and T > 0.5 whereas for 0.25 < T < 0.5 it forms a AB-CDW. 
The phase sequences which occur on increasing temperature can also be very sophisticated. 
In fig. 11(a) one sees that the system at 5 ^ 0.29 undergoes the transitions HOM — yPS- 
Ab/hom — ^AB-CDW — )-HOM on raising T from zero to 1. Another example is found in fig. 
11(b) at 5 ^ 0.165 where the sequence of transitions on raising T from zero to 0.6 is PS- 
AB /hOM ^PS-AB /iP ^IP ^AB-CDW ^IP S-HOM. 

Summarising the results of the section on finite temperatures we conclude that second 
order transitions become more important at T > 0. Phase separation is shifted to higher 
values of the interaction. The AB-CDW is stabilised for finite doping and exists not only at 5 = 
0. Equally the relevance of the incommensurate phase is enhanced. For U G [C/jph, t^ip/ps] = 
[1.914,4.212] and for U G [0, C/jpl] = [0,0.572] the incommensurate phase occurs at finite 
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temperatures but not at zero temperature. Additionally we found re-entrant behaviour as 
function of the temperature in the sequence of phases at given doping and interaction. 

5 Discussion 

In this paper we discussed the phase diagram of the model of spinless fermions with repulsive 
nearest neighbour interaction on a hypercubic lattice in infinite dimensions. In the limit 
d^oc the perturbation series can be resummed and one finds that the Hartree approximation 
becomes exact. This allows investigations in great detail of the phases and their spatial 
structure. In the whole parameter space of interaction, doping and temperature we examined 
the symmetry unbroken homogeneous phase (hom), the AB charge density wave (ab-CDw) 
and the incommensurate phase (ip). We performed a complete stability analysis for the HOM 
and the ab-CDW relying on the rpa which is exact in the limit d ^ oc. For the spatial 
structure of the incommensurate phase we proposed a very general ansatz (23) in product 
form. We were able to show that within this framework the free energy of the IP depends 
only on the amplitude of the order parameter and on the parameter of incommensurability 
7] which is directly linked hy rj = 2h — 1 to the relative frequency h that the order parameter 
has the same sign on two randomly chosen adjacent sites (see appendix A). 

At zero temperature we found that the system is in the HOM phase for large dopings 
and enters the AB-CDW on decreasing doping. For small values and for large values of the 
interaction this transition is of first order and implies a region of phase separation between 
the HOM phase and the AB-CDW (pS-Ab/hom). In the intermediate range U E [f^iPL, t^ipn] = 
[0.572, 1.914] the sequence of transitions for 5 — )■ is more sophisticated since the system 
enters first the IP via a second order transition (see fig. 6). Then the system undergoes a first 
order transition iPi — )-ab-CDW with phase separation between the IP phase and the AB-CDW 
(ps-ab/ip). 

Finite temperatures favour second order transitions; the first order transitions and the 
regions of phase separation are shifted to larger values of the interaction. For T > the 
ab-CDW is stable for finite dopings, too. A very important result is that the region of the 
appearance of the IP is enlarged by finite temperatures. Besides the two important interaction 
values [/jpL and [/jph a third one, C/jp/ps = 4.212, must be taken into account. For values of 
the interaction below f/ip/ps the IP appears in a certain temperature range. If U is not in 
the interval [C/ipl, f^ipn] this temperature range does not extend down to T = 0. The three 
important interaction values C/jpl < C/iph < t/ip/ps divide the positive U -scale into four intervals 
of different generic behaviour. The generic phase diagrams are given in fig. 11. 

Our results can also be used to gain a better understanding of the phase diagram of other 
models such as the Hubbard model and the extended Hubbard model. The argument why this 
transfer is reasonable comprises two points. The first is the similarity of the weak coupling 
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expansion (Hartree approximation) of the (extended) Hubbard model to the treatment of 
spinless fermions presented here. Of course, there are even more phases present in the Hubbard 
model since there are charge and spin degrees of freedom. But the weak coupling calculation 
for charge density waves or spin density waves of a given wave vector are identical to the 
calculations for the corresponding charge density wave in the spinless fermion model except 
for a factor of 2 due to the presence of two fermion species in the Hubbard model. The second 
point concerns the influence of quantum fluctuations. They are suppressed in the spinless 
fermion model as treated here by the limit d— )-oo but they are present in the Hubbard model. 
Fortunately van Dongen showed that the effect of these fluctuations can be accounted for 
in the weak coupling regime by a renormalisation q of the order of 1 [52, 53]. Moreover, he 
showed that in the vicinity of half filling this renormalisation is the same for the different 
phases [45]. From these considerations we conclude that the Hubbard model with pure on- 
site interaction and the Hubbard model with pure nearest-neighbour interaction display in 
the weak coupling limit the same phase diagrams as we calculated except for a rescaling of 
the doping and of the interaction. The general extended Hubbard model should have a phase 
diagram which is in many aspects similar to the ones shown in the present work. 

For instance, van Dongen's result 0.56 = C/„„_aite + 2C/„e„^st_„^igi,bour in our units, i.e. the 
interaction values above which the IP is more favourable than the PS-Ab/hom, is in perfect 
agreement with C/jpl = 0.572 up to the factor 2 due to the spin degeneracy. The difference 
between 0.572 and 0.56 results from the fact that van Dongen compared the analogues of the 
asymptotic expressions (13) and (22) whereas we define [/jpL from the numeric results in fig. 
6. 

At small but finite interactions the renormalisation due to quantum fluctuations is not only 
a simple constant. Nevertheless the qualitative picture remains the same as can be seen in a 
comparison to QMC results by Freericks and Jarrell [50]. They define the temperature ratio 
Tj/Tjy of the temperature Tj at finite doping where the IP sets in and the Neel temperature 
Tjy at half filling where the AB phase (here: spin density wave) sets in. This ratio is found 
to be fairly constant at about 0.57. It is in good agreement to our results for the ratio of 
the temperature co-ordinate of the points a in figs. ll(b)-(d) to the temperature where the 
AB-CDW appears at half filling. We find 0.31 at U = 2.5, 0.51 at [/ = 1 and 0.56 at U = 0.45 
(from figs. ll(b)-(d)). Naturally the agreement becomes excellent for smaller values of the 
interaction. In the comparison of interaction values a renormalisation due to fluctuations of 
the order of 3 should be kept in mind [50]. 

In contrast to the results of van Dongen, Freericks and Jarrell find incommensurate order 
and no sign of phase separation (although they cannot rule out the latter completely). These 
seemingly differing findings can be reconciled by our results. In the quantum Monte Carlo 
calculations the phases at T = are found by extrapolation from finite temperature results, i.e. 
one investigates which phase occurs first on lowering the temperature from the paramagnetic 
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phase (corresponding to the HOM phase). Looking at figs. ll(b)-(d) one notes that it is either 
the AB phase or the incommensurate phase that is found. There is (almost) no chance to find 
a sign of phase separation, especially since the divergence of the compressibility occurs behind 
the right phase separation line. Adopting, however, van Dongen's view point at T = and 
U ^ fig. 6 clearly shows that AB phase separation (PS-AB /hom) dominates. This resolves 
the seeming contradiction and adds to the understanding of the phase diagram in the more 
complicated case of the Hubbard model. The QMC result confirms our statement that finite 
temperatures favour the appearance of incommensurate phases also in the Hubbard model. 
We are convinced that the infinite dimensional Hubbard model displays also phase separation 
between the antiferromagnetic phase and an incommensurate spiral phase corresponding to 
PS-AB /iP. 

Further investigations should incorporate the effect of fiuctuations which implies the con- 
sideration of 1/d corrections for the spinless fermion model. This was done so far only at half 
filling [48] due to the high complexity of the phase diagram at finite doping. 

We thank D. VoUhardt and P. van Dongen for helpful discussions. We are indepted to 
V. Janis for valuable advice. Furthermore we gratefully acknowledge the financial support 
of the Deutsche Forschungsgemeinschaft in the Sonderforschungsbereich 341 (GSU) and the 
Stichting Fundamenteel Onderzoek der Materie (RV). 



Appendix A 

In order to derive that the ground state energy per site, more generally the free energy F 
per site, depends only on and h if the order parameter b(r) is assumed to be given by the 
ansatz (23) we consider the expansion of F in powers of b(f). The direct interaction part 
(representing TTrSG — $; cf. first two terms in (30)) is given by 

= (C//2)[(l-/i)(n + 6o)(n-6o) + (V2)((n + M' + («-M')] (34a) 
= (C//2) (n^ - (1 - 2h)bl) = (C//2) (n^ + r]bl) . (34b) 

The kinetic part (representing — Trln(— G); cf. last term in (30)) results from an expansion 
in the inhomogeneous chemical potential given by Ub{f) 

oo Tjn 

= TJ2Y.—Sn with 

n-l 

^» ^ 5] K'^i)9{ri-r2)b{f2)g{r2-f-^)...b{fjg{f^-f^) (35) 

where g(f) is the propagator in real space. Its frequency dependence is omitted since all 
propagators are used at a fixed frequency due to the static character of the Hartree 
theory. Based on the following representation of the propagator g{uj^, k) in momentum space 

9{<^n^k) = {iu^ - ej: + ft)-^ = -i dXexp - ej: + jl)Xj (36) 
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for o;^ > (for o;^ < one uses iX — > —iX) the propagator g(f) in real space can be written 
as 



9{r) 



g{u^, k) exp [ikrj = -i dX exp (-Aw„ + ifiX) JJ /(r*^\ A) 



27rY 
^ dk 



/(r, A) := / — exp (iX\j2ldcos[k) + ikr\ 



= 5{r) + i{X/2)J2/diSir + 1) + 5{r - 1)) + 0{d-^) . 



(37a) 

(37b) 
(37c) 



The superscript in r^^^ stands for the j-component of the vector r. The representation (37) 
together with (23) allows to write the essential part of as a product 



5. = bl 



dX^... dA„e-<"»-^'^^<^i+-+^»'(-i)"'P("'({A.}) 



(38a) 
(38b) 



L 



I ^ ,/ 2 vi' ji 



where the system size iV = L*^ is split up assuming a hypercubic shape of the system. The 
expression in (38) is still exact in any dimension. Inserting the expansion (37c) one finds for 
odd n that already the leading order d^ vanishes since b(f) vanishes on averaging over the 
lattice {l/N)J2fb(f) = 0. For even n the leading order of (37c) yields 1. Now we apply an 
argument for d ^ oc similar to the one used in the derivation of the central limit theorem 
expanding the logarithm of P^^\{X-}) 

ln(P*"H{Aj))=Eln^-^ E E \X,v^-''''{r^u^;^{r,+ a) + 0{dA (39a) 

i=l \ ■^'^^ r^;a^±l n>I'>l>l J 

= - E \\>{h + il-h){-lf-^^)+0{d-'). (39b) 

n>l'>l>l 

In the last step we used = \ for even n and u^'~^' {r^)u^: ~\r-^ + a) = \ for even I' — I. For 
odd /' — I the expression (r^)'U^^^(r^ + a) takes the value 1 with the relative frequency h 
and the value —1 with the relative frequency 1 — h. 

The crucial point in (39b) is that the result depends for d ^ oc only on h and not on 
further details of the functions Uj. Via (38) the free energy functional depends only on b^ and 
on h (or equivalently on rj = 2h — 1). 

We turn now to the explicit derivation of (19). Starting from 

-1 



G{iuj^,k]h) = 



-A 



(40) 



with A := (1 — 2/i)6q[/ = —rjbQU, it is appropriate to split = '^^ + '^^ with 



hd 



(41) 



n = l 



n = hd-\-l 
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so that e- - = — holds. The one-particle energy resultmg from (40) is e(u,v) : = 



u + sgn{v)^v'^ + A^. The distribution of u and v is again easily found via the central limit 
theorem which permits to compute the DOS pfj{uj) 

. /^ii^.,(_[.V...V(i-.)]/2).(.-.(.,.,). ,«) 

-co 27r^/i(l - h) 

Performing the -u-integration and substituting h by (r] + l)/2 yields an expression for p^{uj) 
which leads together with (34b) to the free energy 



F(n,A;r?) = - -n^ + i^n + T J du \n {1 - Mco)) p^{lo) . (43) 

— oo 

Note that f-p{uj) = 1/ [1 + exp(/3(a; — jl))] In the limit T — )■ one obtains the ground state 
energy in (19). 



Appendix B 

The matrix elements in (29) are calculated from a single dressed bubble using (17) 

>2 



where we write iu^ for ^ti;^ + fl for the sake of shortness. The sums run over half the 
Brillouin zone (e- < 0) in order to avoid double counting. Performing the sum of Matsubara- 
frequencies uj^ and using (4) for the /c-sum yields 



iv) = vp / de^de2p{e^,e2;'n)^^—^p^ — 2^ — (/f(-z') - /f(^')) (45a) 
(^) = vp / de^de^p{e^,e^\r]) ^ 2 ' (45b) 

./ 69 — €1 



where vp stands for the principal value and z' := ye^ + A^. The partial fraction expansion 

^pf^^J = vp(-)vp(— ^ + ^^) -7r25(ei) (5(61-62) + 5(61 + 62)) (46) 

is useful to simplify the integration in (45). After the substitution 62 = u^2{l — rp') + rje-^ the 
ti-integration is mapped on the representation (26) of the Dawson-function and one obtains 
finally 

Xtiv) = 1 "^/7'^1 [.^D(a._) + A^D(^,)]F_(.)d^--^F_(A) (47a) 
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A 



/ 



ttA 



B{u_) + Biuj^)\F4z)duj-^-j===F4A) with (47b) 



a;^ := uJ{l±v)/{2{lTv)) F ^{z) := M- z) ± Mz) z := Vu^ + . (47c) 
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Figure captions 

Figure 1: First three contributions to the Baym/KadanofF generating functional 
$: (a) leading order d^; (b) and (c) order d~^. 

Figure 2: Diagrammatic representation of (a) the Dyson equation (b) the general 
Bethe-Salpether equation and (c) the rpa series for the density-density response 
function x(A;) resulting from (b). 

Figure 3: (a) solid curve: interaction values where x(?y„p,) diverges at T = 0; dotted 
curve: interaction values where x(— 1) diverges, (b) the corresponding values of 
r?„pt, for < 5 < (5end = 0.308 holds r/,^, > -1. 

Figure 4: Ground state energies at U = 1.2, for the sake of clarity a linear func- 
tion of S is added (a « 0.5561, « 0.1946); dotted curve : homogeneous phase 
(hom); dashed-dotted curve: AB charge density wave (ab-CDW); solid curve: in- 
commensurate phase with optimised rj = rj^^t (iP)- Divergence of x(— 1) at ^jj^.^B! 
divergence of x(^opt) S^-^_jp. The Maxwell construction for the phase separation 
HOM-^ — ^AB-CDW is a straight line between the points at S = 5ps_ab and 5 = 0, 
the Maxwell construction for the phase separation IPi — ^AB-CDW is a straight line 
between the points at 5 = 5ps_ip and 5 = 0. The pure IP is present between 5ps_ip 
and 5di^_ip. 

Figure 5: Phase separation between the HOM phase and the ab-CDW at T = 0. The 
phase separated region is marked ps-ab/hom. Dotted curve: interaction values 
where x(~l) diverges; solid curve: PS-AB/hom boundary. The pure AB-CDW is 
only stable at half- filling (5 = 0). 
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Figure 6: Phase diagram at T = 0: HOM: homogeneous phase, AB-CDW: AB 
charge density wave, IP: incommensurate phase, PS-AB/hom: phase separa- 
tion between HOM and AB-CDW, PS-AB/ip: phase separation between IP and 
AB-CDW. Sohd curves: first order transition HOM-^ — )^AB-CDW; dashed curve: 
first order transition IP-^ — )-AB-CDW; dashed-dotted curve: second order transi- 
tion HOM^ — S-IP. 

Figure 7: (a) sohd curve: interaction vahies where xiVopt) diverges at T = 0.1; 
dotted curve: interaction values where x(— 1) diverges, (b) the corresponding 
values of r?„pt, for S G [5start,4nd] = [0.055,0.319] holds r]„p, > -1. 

Figure 8: Relevance of the incommensurate phase: the dashed and the dotted 
curve enclose the region where rj^^^ > — 1 is possible. The dotted curve shows 
the points of the continuous transition r]„pt — )■ — 1 on leaving the enclosed region 
whereas the dashed curve depicts the points where the transition 7]„J,^ —1 is 
discontinuous. The solid curve depicts the S and T co-ordinates of the critical 
points of the AB phase separation. The co-ordinates of the intersection a define 
'^ip/ps cind Tjpypg. 

Figure 9: Generic phase diagram (solid lines) at finite temperature T = 0.1 < 
Tip/pg. Besides the nonsegregated phases AB-CDW, IP and HOM regions of phase 
separation occur between the AB-CDW and the incommensurate phase (ps-ab/ip) 
and between the AB-CDW and the homogeneous phase (ps-ab/hom). The point 
b is tricritical. The ps-ab/ip starts there: curves above point b mark first order 
transitions. The curves ab and ac mark second order transitions. 
The dashed and the dotted line show the phase diagram if the new ansatz for the 
IP were not taken into account. The dotted curve marks a second order transition 
from AB-CDW to IP calculated from (28). Note the significant difference between 
this curve and the curve ab. The dashed curve marks the beginning of a first 
order transition from AB-CDW to HOM. 

Figure 10: Generic phase diagram at finite temperature T = 0.2 > Tip/pg. Solid 
curve below the critical point a: second order transition HOM^ — ^AB-CDW; solid 
curves above point a: first order transition HOM-^ — )-AB-CDW, left branch U{S-^), 
right branch U{S2) (for 8-^^,62 see explanation to eq. (32)). 

Figure 11: Generic phase diag rams at fixed interaction, (a) U — 5.0 > f/jpypg: 
point a is the critical point where the second order transition line (above a) and 
the first order transition lines (below b) meet. 

(b) [/ = 2.5 G [C/,ph,C/:p/ps]; (c) U = 1.0e [C/ip„C/iph]; (d) C/ = 2.5 G [0,U,^^]. 
The points and the lines between them in figs, (b)-(d) have the same meaning as 
in fig. 9. 
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